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Summary 


How to calculate the visible length and ground-level concentrations resulting from a wet polluted plume 
emanating from an elevated stack is illustrated. In this example the stack is assumed to be on the outskirts of a 
fictitious town called Midville. It operates only during daytime, ona small slope within a housing area — typical of 


many such installations. 


1. Introduction 

It is assumed that an authority in Midville (somewhere 
in the Midlands of England) is proposing to build a new 
municipal incinerator with a stack for waste gases. The 
gases will contain sufficient water vapour to form a 
visible condensed plume which will be of potential 
concern to the local inhabitants. 

The aims of the analysis in this paper can be 
summarized as: 

(a) To determine the visible length of the plume in 

different weather conditions. The higher the humidity 

of the ambient air and the lower its temperature, the 

more slowly do the fine drops in the stack plume 

evaporate. Wind speed is also important since the 

stronger the wind the sooner, in terms of downwind 

distance, does the plume evaporate. 

(b) To determine the maximum ground-level con- 

centration of each significant pollutant contained in 

the plume, as measured over several minutes — say 

10 minutes, wherever it may occur, as a function of 

stack height and weather conditions. 

(c) To determine where these maximum concen- 

ations occur. 

(d) To determine whether or not any surrounding 

houses are at risk from the effluent gases and what 

remedial action needs to be taken. 


2. A list of the basic parameters 
Values of the key parameters will be given the 
following values by way of illustration: 
The temperature of the effluent gases at the stack 
top =140 °C. . 
The efflux velocity (determined by afan)=15ms_'. 
Maximum gas flow rate = 600 N mh". 
Internal stack diameter = 370 mm. 
Water mixing ratio at the stack top = 0.3. 


The stack height can be chosen to give least nuisance, 
remembering that local inhabitants will tend to object to 
a tall unsightly stack even more than to the occasional 
odour. Calculations are made for 16 m, 20 m and 23 m. 
Extreme values of 10 m and 30 m are also considered. 

The effluent gas will contain a mixture of substances 
including hydrochloric acid (HCl), sulphur dioxide 
(SO2), nitrogen oxides (NO,) and carbon monoxide 
(CO). The assumed emission rates are: HCl 300 gh", 
SO2/NO,/CO 600 gh". 


.3. The visible length of the plume 


The gases rising up the stack contain a considerable 
amount of water which rapidly condenses on leaving the 
stack, giving a visible plume. As the plume broadens by 
turbulent mixing with the environmental air, drier and 
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colder air is mixed with the moist, stack air. This cools 
the plume, tending to retard its evaporation, but also 
dries it. How long the plume remains visible before this 
drying-out process results in total evaporation depends 
on the temperature and humidity of the environmental 
air. It also depends on the amount of natural turbulence 
in the air, which in turn depends on the stability of the 
lower atmosphere. The more stably stratified the air 
becomes, the lower is the level of turbulence and the 
slower the drying out of the plume. 

The equations governing this process are now set out: 

Let a= the emission of air = 1.67 X 0.00127 X 10° gs' 

= 2.12 X 10° gs", and 

b = the emission of water = 0.3 X 2.12 X 10° = 0.64 X 

10°gs'. 

Define s = moisture content of the environmental air 

as a mixing ratio ingg. 


The total water emission is thus (b + as) g s ', and the 
total dry air emission is (I—s)a gs‘. 

The specific heat of air at constant pressure = 0.24, 
and that of water vapour = 0.48. 

If, at a distance x downwind of the stack, the mass of 
entrained environmental air = m g, if T;= the exit gas 
temperature = 140°C, 7,=the temperature of the 
environmental air, 7= the temperature of the plume at 
x, and $= the mixing ratio of the plume at x, then: 

H;= heat content of the flue gases = 

(0.24(1 — s)a + 0.48(6 + sa))T;, where a = 2120, b = 640, 

and 7;= 140. 

W,;= the water content of the flue gases = (b + sa). 

H.=the heat content of the entrained air = 

(0.24(1 — s)m + 0.48 sm)T;. 

W.=the water content of the entrained air=sm 

(which is a function of s and x). 

H, = the heat content of the plume = 

T (0.24(1 — s)(a + m) + 0.48(6 + sa + sm)). 
W, = the water content of the plume = (a+ b+ m)S. 


Now H,= H;+ H., so that ((1 + s)a + 2b)(T;— T) = 
(} + spel TF — #). 

Similarly W,= W;+ W., so that (a+ b)(S;— S)= 
m(S—s), where S; is the water content of the flue 
gas = 0.3. 

If we eliminate m from these two equations we arrive 
at the ‘Mixing Curve’ equation: 


(1+s)a+2b (B= r) - (22) 
stb (Ss) “8 §-,)- 


Now after some mixing has occurred, certain simplifi- 
cations can be made to this equation. We may assume 
that T7;>T, S->S and s<l. The Mixing Curve 
equation then simplifies to: 








S—s=A(T— T,) 


where A = (a + b)S;/(a + 2b) T;= 0.00174 if the values 
previously given are inserted. Strictly this value of A 





applies when T= 0. At T= 10 a more accurate value of 
A is 0.00187, and at T= 20 A is 0.00203. Now s is the 
mixing ratio of the environmental air and may be 
expressed as rs-(7,), where s- is the saturated mixing 
ratio at temperature 7,. If T, is expressed in degrees 
Celsius then a good fit to s.(T,) is given by: 


s{T) = 10°°(3.839 + 0.3008 7+ (0.005567 + 0.0003587T") 


r is the relative humidity (expressed as a fraction, not as 
a percentage). 

The value of the plume temperature when the plume is 
on the point of evaporation is found by an iterative 
process. To begin, assume T= T,. Insert this into the 
equation for s. and find s.. Insert this into the simplified 
mixing curve equation and find a new T. Repeat again 
and again until the solution converges to a constant 
value of T, and a corresponding value of S. 

To find at what distance this occurs we must define 
how rapidly the plume broadens with distance. During 
the daytime, the plume will only persist a significant 
distance downwind in very humid conditions (when r is 
only a little less than 1). Such conditions occur in neutral 
stability conditions when the instantaneous plume 
width grows like: 


R(x) = R(O) + 0.176x 


where R(x) is the radius (in metres) and x is the 
downwind distance, also in metres. The formula is 
consistent with the plume acting as a cone with a half- 
angle of about 10°. To find R(0), we equate the volume 
flux emerging from the stack with the effective volume 
flux assuming the plume reaches its effective source 
height very rapidly and thereafter behaves like a 
horizontal plume being carried by the wind: 


a+b 
p 





UrR’(0) = 


where p = the air density = 1270 gm’. 
Remembering that m is the mass of entrained air, an 
earlier equation gives 


_ (at bXSr— S) _ 
vali ps “ 


0.3(a + b) 
0.0019(T— T:) 


approximately. The total mass of the plume M= m+ 
(a + b) = 12707 R’U. These equations relate the plume 
temperature T to distance x, and hence tells us the 
distance when the plume evaporates. These distances are 
given in Table I. 

The figures summarize the results. Fig. | shows the 
saturation curve as a function of air temperature and its 
mixing ratio. The graph also shows a sample of ‘mixing 
curve’ lines which are virtually straight over the range of 
temperatures shown. Thus a plume from the stack 
emerging into an airstream whose temperature is 10 °C 
and relative humidity is 0.5 (as a fraction), evaporates 
when its temperature has fallen to 13 °C. 
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Table l. A table of the distances x, in metres, at which the plume from the 
stack evaporates, as a function of ambient temperature 7,, relative humidity r, 
and wind speed U (m s ') measured at 10 metres above ground 

















T, r Tevap U=1 U=3 U=5 U=7 U=10 
0 0.4 1.45 44.9 25.9 20.1 17.0 14.2 
0.6 0.96 56.0 32.3 25.0 21.2 17.7 

0.8 0.48 81.1 46.8 36.3 30.6 25.6 

0.9 0.24 116.6 67.3 52.1 44.1 36.9 

0.95 0.12 166.8 96.3 74.6 63.0 32.7 

0.98 0.05 266.4 153.8 119.1 100.7 84.2 

5 0.4 122 35.4 20.4 15.8 13.4 11.2 
0.6 6.47 44.4 25.7 19.9 16.8 14.1 

0.8 Ye | 64.8 37.4 29.0 24.5 20.5 

0.9 5.36 93.7 54.1 41.9 35.4 29.6 

0.95 5.20 134.4 77.6 60.1 50.8 42.5 

0.98 5.07 215.3 124.3 96.3 81.4 68.1 

10 0.4 13.53 27.2 15.7 12.2 10.3 8.6 
0.6 12.31 34.6 20.0 15.5 13.1 10.9 

0.8 11.14 51.1 29.5 22.9 19.3 16.2 

0.9 10.57 74.4 42.9 33.3 28.1 235 

0.95 10.28 107.2 61.9 48.0 40.5 33.9 

0.98 10.11 172.4 99.5 77.1 65.2 54.5 
15 0.4 21.22 19.5 11.3 8.7 7.4 6.2 
0.6 18.93 25.6 14.8 11.4 9.7 8.1 

0.8 16.88 38.8 22.4 17.4 14.7 12.3 

0.9 15.92 57.3 33.1 25.6 21.6 18.1 

0.95 15.46 83.2 48.0 37.2 31.4 26.3 

0.98 15.18 134.5 77.7 60.2 50.8 42.5 

20 0.8 23.52 a3 15.7 12.2 10.3 8.6 
0.9 21.66 41.7 24.1 18.6 15.7 13.2 

0.95 20.81 61.7 35.6 27.6 23.3 19.5 

0.98 20.32 100.6 58.1 45.0 38.0 31.8 

25 0.9 28.90 25.8 14.9 11.6 9.8 8.2 
0.95 26.70 41.3 23.8 18.5 15.6 13.1 

0.98 25.65 69.1 39.9 30.9 26.1 18.5 

0.018 
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Figure 1. The saturation mixing ratio is a function of air temperature. Below the curve, the air is unsaturated. The nearly-straight lines 
running up to the saturation curve are the mixing curves for different ambient air conditions, but the same emerging plume conditions. 
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Fig. 2 shows how the maximum downwind extent of 
the visible plume varies with ambient air temperature, 
relative humidity and wind speed. It assumes near 
neutral stability conditions, which typically occur about 
70% of the time in the United Kingdom. In unstable 
convective conditions, the plume would evaporate 
earlier whereas in stable conditions it would take longer. 
However, stable conditions are likely to occur during 
the daytime mainly at the very beginning and the very 
end of the working day when the stack is not operative. 
If stable conditions occur at other times (perhaps when 
milder air replaces bitterly cold air in winter) then the 
plume could be a temporary nuisance, and a case could 
be made for suspending incineration during such a 
period. During these times the plume length could be 
increased fivefold above the values given, although this 
is somewhat speculative. 





As will be seen from the data, the visible length of the 
plume in neutral conditions exceeds 100 m in length 
only in very high humidity conditions and at low wind 
speeds and temperatures. 


4. Maximum ground-level concentrations 
Since concentrations are directly proportional to 
emissions, the calculations will be made using an 
emission rate @= 1 gs '. The ground-level concentrations 
are, however, not only dependent on source strength but 
on the effective height of the source, the wind speed and 
the stability of the air. The effective source height is 
defined as that height which, if the influences of gas 
temperature and efflux velocity were to be ignored, 
would give the same ground-level concentration pattern. 
The efflux gases are emitted within a vertical jet 
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Figure 2. The maximum downwind extent of the visible plume from the stack, with different wind speeds and different ambient relative 


humidities and temperatures. 
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travelling upwards at 15 ms‘ and at a temperature of 
about 140 °C. Both these factors make the plume rise to 
some effective height A. significantly greater than the 
stack height A. 

The plume rise Ah = h. — A. can be calculated using 
formulae given by Briggs (1984). Briggs states that the 
plume rise is given by: 


Ah= (“43 + 3fex. ) P 
BU 2B U 
where U= wind speed inms', W= 15 ms‘, the efflux 
velocity, x = the downwind distance in metres, Fi» is the 
momentum flux, F; is the buoyancy flux, and B is given 
by (0.44 1.2U/ W. 

However, Briggs suggests replacing B in the second 
term (the buoyancy term) on the right-hand side by 0.6 
to be consistent with the asymptotic behaviour of 
buoyant plumes. 





The first term in the brackets on the right-hand side of 
the equation is given by 


(3) (eatéare) 7 
4) \@2+06/R)/ ~*~ 
where R= W/U and W=15ms'". 

The second term is given by putting 


P ie ack ? 
Fo =e — Wr 
where W = 15, r, is the stack internal radius = 0.185 m, 


Tstack = 140 + 273, Tair = 15 + 273, g=9.81. 

If h; = 23 m, then A; is given in Fig. 3 as a function of 
wind speed U. The effective height h. increases rapidly as 
U falls below about 3 ms '. Increasing or decreasing the 
stack height simply raises or lowers the effective height 
by the same amount. 


100 T T T T T T T T T 
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Figure 3. The effective height of the plume for a 23 m stack when 
the combined effects of the high temperature and the efflux velocity of 
the emerging plume are taken into account. Otherwise, effective height 
is principally a function of wind speed. 


The plume rise may be ultimately halted by an 
elevated stable layer through which the plume cannot 
penetrate. However, before that, the plume may 
effectively level out and diffuse very much like a passive 
plume. It is assumed that, for modest emission rates, this 
occurs at 100 m downwind from the source. Since this is 
likely to provide an underestimate of the true plume rise, 
it will overestimate the ground-level concentrations to 
some extent. 

Having so defined h., the plume is then assumed to act 
like a neutrally buoyant plume emanating from h-.. The 
approach of Weil (1988) is used. It is a very simple and 
effective method for neutral and unstable conditions. 
Essentially it calculates the concentration by assuming 
the plume elements travel along straight lines radiating 
out from the effective source according to the probability 
distribution of the turbulent velocities at the effective 
source. Elements are assumed to reflect off the ground 
or any elevated capping inversion above the source 
rather like light rays would off reflecting mirrors. 
Although this is not a perfect representation of what 
happens in the atmosphere, field trials indicate that the 
model gives reasonably satisfactory estimates of measured 
concentrations. Normally, out to the distance where the 
maximum concentration occurs, the effect of any 
capping inversion is not felt at ground level (except 
indirectly in its subtle effect on the structure of the 
turbulence). 

The plume spreads not only vertically, of course, but 
horizontally as well. The acrosswind spread is particularly 
important in lowering the average concentration of 
pollutants in the plume with downwind distance. 
Following Briggs (1985) we assume the acrosswind 
distribution has a Gaussian shape: 





Cy) = exp( — a) 
0; 


V2, 
where o, (the so-called ‘standard deviation’) is a measure 
of the width of the plume and is approximately given by: 


oy_ _0.6X 

h- (1+2X)" 

where X = wsx/(Uh), we is the convective velocity = 
(gHh/pcpT)"3, and H is the surface sensible heat flux in 
W m’”, his the depth of the mixing layer. 

Typical values of h, H and ws can be ascribed to pairs 
of wind speed U and stability category P (see Table II). 
The mixing depth A strictly depends on a number of 
factors which are difficult to assess in any given situation 
without detailed meteorological data both at the time of 
interest and back in time along the trajectory of the air 
mass. Since such information is unknown in the present 
analysis, a rather simple approach which may introduce 
small errors into the results has to be used. It is believed 
that these errors do not significantly alter the conclusion. 
This belief is based on a sensitivity study by which 
changes in the results are obtained by changing the 
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Table Il. 


A table of parameters governing the ground-level concentrations downwind of 


the stack. H is the heat flux in W m °, wg is the convective velocity in m s |, his the estimated 
mixing depth in metres and uy is the friction velocity in ms '. The plume rise is given in 
metres as is the effective height of the source (= plume rise + stack height (23 m)). 





Wind speed 1-3 kn 4-6 kn 7-10 kn 11-16 kn 2 17kn 
(average) Ims' 25ms' 4.25ms' 6.75ms' >9ms' 
Plume rise (m) 62 26 18 14 13 
Effective h (m) 85 49 41 37 36 
P Param. 
A H 150 300 
We 1.86 2.63 
h 1634 2368 
Ux 0.1476 0.2928 
B H 20 100 250 
Wy 0.680 1.52 2.40 
h 651 1459 2316 
Ux 0.1176 0.2647 0.4272 
£ H 10 40 110 260 
Wy 0.481 0.961 1.594 2.451 
h 495 1028 1699 2618 
Ux 0.110 0.2471 0.4056 0.6275 
D H 0 0 0 0 0 
Wy 0 0 0 0 0 
h 217 543 923 1466 1955 
Ux 0.0869 0.2172 0.3692 0.5863 0.7821 





values of h within sensible limits around the values given 
by our simple method. 

‘The simple method uses the formula for the develop- 
ment in time of the depth of the convective mixing layer 
by Carson (1973): 


dh’ _ 2(1+2A)H@, 1) 
dt pcpy(t) 


where H(0, £) is the surface sensible heat flux at z= 0 at 
time ¢, A is a factor to allow for the entrainment of heat 
from above the mixing layer and is usually given a value 
of 0.25, y(t) is the potential temperature gradient in the 
air above the capping inversion, and is given a default 
value of 0.003 °C m'. The equation can be solved for h if 
an analytic form for HA(0, 2) is given. In steady 
meteorological conditions H(0, t)= Hnsinwt where Hn 
is the peak midday value of the heat flux, and is the 
inverse time-scale of the diurnal cycle. The maximum 
value of / is obtained in the afternoon, and it is this value 
that is used. It turns out that A depends almost entirely 
on H,, in a very simple way: h = 130H nm. In the morning 
the depth of the mixing layer is likely to be less which 
leads to weaker turbulence and lower ground-level 
concentrations until the physical limits to vertical 
spreading imposed by the lower depth push up the 
concentrations to some degree at distances beyond 
which we are concerned. 





So as to obtain the optimum estimate of h, and also to 
blend in with alternative estimates of A in neutral 
conditions, the value of h=hc given by Carson’s 
equation is combined with h = (ux/4f) commonly used 
in neutral conditions: 


hr = het (us/4f/. 


The values of H are obtained from Smith’s version of 
the Pasquill scheme for estimating stability. Fig. 4 
shows Smith’s curves which relate P, Uand H. The table 
of parameters (Table II) also gives the convective 
velocity ws and the friction velocity ux. The convective 
velocity has already been defined and is used in the 
non-dimensional form of the downwind distance. The 
friction velocity, which is used in the mixing layer depth 
in neutral conditions, can be inferred from the 10m 
wind speed and the heat flux by the following equations 
(which are standard equations used in micrometeorology 
and won't be explained properly here): 


- O&— DO. +t 1) 
Ux = 0.4U/ finde s+ 


+ 0.0349(tan 'y — tan’! vo) 


where y=(1+ 160/L)“, yo=(1 4+ 1.6/L)%, L is the 
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Figure 4. The revised scheme for Pasquill Stability based on heat flux and wind speed. 


Monin—Obukhov length-scale = 94 700u%/ H, and the 
surface roughness has implicitly been assumed to be 
0.1 m. In neutral conditions us = 0.087U. To solve for 
Ux in convective conditions the neutral value of Ux is 
found, this is inserted into the expression for Lto find L 
(for given heat flux H), y and y. are found and inserted 
into the convective expression for us. This new value of 
ux is then reinserted into the L expression and so on 
until a converged value of ux is reached. To determine 
the ground-level concentrations, the formula provided 
by Weil (1988) will be used. The concentration below the 
plume centreline is given by: 





7 a(X) { _ (a@/ X — 0.35)" 
C(X) & y 1.841 pexp( ais dt 


2 
+ 0.361 exp(— +" As ) 


where a(X) = 0.665 ((1 + 2X)*/X). 
Q is the emission rate in gs', a=h,/h, and X= 
(Wx/ Uh). 


5. Results 

The equation for C(X) above has been solved for 
Q=1gs' and for all the combinations of U and Pin 
Table II, and for five stack heights: 10 m, 16 m, 20 m, 
23 m and 30 m. 

In each case the concentration distribution has the 
following characteristics: on the ground the concentration 
is zero close to the base of the stack but after some 
distance downwind it suddenly rises quite rapidly to its 
peak value, after which there is a gradually decreasing 
tail which persists for a long distance downwind. As an 
example, the results for stability P= Band a wind speed 
of 2.5ms_‘ are given: 


X=0.01: C=0; X=0.02: C=0; X= 0.03 : C= 3.9; 
X= 0.04: C= 27.2; X=0.05: C= 47.7; X=0.06: 
C=52.1; X=0.07: C=47.6; X=0.08: C= 40.5, 
etc. 


Two important conclusions can be drawn from the 
calculations: 

(i) C(X) is at a maximum with wind speeds around 
2.5ms'. For example with P= C the maximum 
concentration varies with wind speed as follows: 
U=1:C=52.0; U=2.5: C=53.1; U= 42.9: 
C= 42.9; U= 6.75: C= 29.3; 

(ii) C(X) varies only very gradually with stability. 
Cmax is the overall maximum concentration, regardless 
of where it occurs. It is given in Fig. 5 as a function of 
the most important parameter, h,, the stack height. 
Note that decreasing h, from 23 m to 16 m increases 
Cmax from 59 to 79 ug mm’ with Q= 1 gs‘. (With the 
proposed stack, having an emission rate of only 
0.166 gs for sulphur dioxide, the concentration 
increases from 10 to 13 ug m~.) 


Fig. 6 shows the downwind distance xm in metres 
where the maximum concentration is expected. Xm is 
virtually independent of wind speed U, but depends very 
much on stability P and A. 


6. Maximum concentrations around the 
stack 


Ignoring the sloping terrain in the vicinity of the 
stack, the concentrations of some of the pollutants 
emitted from the stack can be estimated. 


6.1 HCl 


The emission rate of HCI is given as 300gh" or 
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Figure 5. The 10-minute-mean peak concentration at ground level 


as a function of stack height for the resulting plume, but with an 
emission rate standardized to 1 gs‘. The units of concentration are 


ugm. 


0.083 gs '. The 10-minute-average maximum con- 
centration, as a function of h,, in ug m ° is: 

h{m) 10.0 16.0 20.0 23.0 30.0. 

Co 6) 66 55 49 33. 


6.2 SO2, NOx, CO 

The emission rate for each of these pollutants is 
estimated to be 600gh' or 0.166gs'. The concen- 
trations are: 

h{m) 10.0 16.0 20.0 23.0 30.0. 

Cmax: 17.4 13.2 11.0 9.8 7.6. 


6.3 The effect of the slope 

It is assumed that the ground rises some 5 metres to 
the north-east, say; and further that a tall house exists on 
the highest ground and is some 8-10 metres tall, so that 
its upper rooms may be some 14 metres above the base 
of the stack. 

The slope of the ground will only be important in the 
relatively few cases of daytime stable conditions when, 
due to the increase of air temperature with height, air 
motions tend to be much more horizontal than in 
neutral or unstable conditions when the air readily 
follows the contours of the underlying terrain. 

Although stable conditions can and do occur at any 
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Figure 6. The distance downwind on flat terrain where the 
10-minute-mean maximum concentration occurs, as a function of 
stack height and Pasquill Stability. 


time of the day if the meteorological conditions are 
right, they most frequently occur soon after dawn and in 
the late afternoon an hour or so before sunset. They tend 
to develop whenever the winds are light and the skies are 
cloud-free. To avoid the great majority of these cases, it 
would be prudent to recommend, in this assumed 
scenario, that the stack should not be operated with very 
light winds after, or before, a cold night, with a south- 
westerly drift moving the air up the slope towards the 
tall house to the north-east, until at the earliest 9 a.m., or 
until the latest 4 p.m. Under these stable conditions the 
upper rooms of the house would otherwise be subject to 
emissions from a stack that is effectively 14 metres 
lower. Calculations show that for a 23-metre stack this 
apparent reduction doubles the maximum concentration 
likely to be experienced, even though the centreline of 
the plume will still easily clear the top of the house. In 
reality the proposed restriction will impose very little 
extra constraint on the stack’s operation since firing 
would not normally take place at either end of the 
working day. 

More precise criteria could be proposed, and a simple 
warning scheme arranged with the Meteorological 
Office at the nearest forecasting office. 
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7. Frequency data of the meteorological 
parameters 

No meteorological observing station exists in Midville 
itself, and the nearest fully operational station is at 
Elmdon. Although some distance away, it is considered 
that the statistics of the relevant parameters are likely to 
be very similar. Between 08 and 18 UTC the Pasquill 
Stability Categories occur with the percentage frequencies 
given in Tables III and IV. All the stable cases (E, F and 
G) occur either at 08, 16, 17 or 18 UTC. The highest 
percentage frequency is seen to be in the neutral D 
category, and over 92% of hours at Elmdon lie within 
C-D, i.e. in neutral or slightly unstable conditions. The 
more unstable conditions tend to occur around midday 
when the incoming solar radiation is at its highest. 


Fig. 7 shows the wind rose for Elmdon, which is less 
than 20 km to the north-west, as well as for Brize Norton 
(further to the south) and Wyton (more than 100 km to 
the east). A very similar wind rose to that for Elmdon 
should apply to the stack. South-west winds are seen to 
be considerably more frequent than other winds. The 
sector from south to west contains 75% more winds than 
the average in the north-east. 

Between 08 and 18 UTC, the wind speed classes have 
the percentage frequencies as set out in Table V. As can 
be seen, very light winds below 2 ms occur in about 
5-6% of all hours, but many of these would be at night 
when the incinerator would not be operating. 





Table lll. 


Percentage frequencies of Pasquill stability categories at Birmingham 


Airport (Elmdon), Brize Norton and Wyton during daytime (08-17 UTC) 




















P A K c D F F  G 
Elmdon 0.38 4.13 24.72 67.39 2.52 0.57 0.28 
Brize Norton 0.32 3.96 22.76 69.73 2.17 0.67 0.28 
Wyton 0.30 3.81 22.33 67.43 2.49 0.50 0.14 
Table IV. Pasquill Stability Category frequencies by hour at Elmdon 1980-89 
Pasquill Stability Category 
Hour A B Cc D D E F G Total 
(UTC) (day) (night) 
0 0 0 0 0 1068 1523 639 422 3652 
l 0 0 0 0 1051 1518 637 446 3652 
2 0 0 0 0 1066 1490 619 477 3652 
3 0 0 0 0 1078 1492 583 499 3652 
4 0 0 0 400 1003 1329 530 390 3652 
5 0 0 0 1160 837 1045 354 256 3652 
6 0 0 21 1739 720 781 229 162 3652 
7 0 2 297 2050 510 549 151 93 3652 
8 0 53 603 2373 261 259 51 52 3652 
9 4 147 927 2574 0 0 0 0 3652 
10 20 211 1227 2194 0 0 0 0 3652 
11 33 258 1348 2012 0 0 0 0 3651 
12 36 268 1315 2033 0 0 0 0 3652 
13 29 233 1207 2183 0 0 0 0 3652 
14 14 189 975 2473 0 0 0 0 3651 
15 3 111 726 2812 0 0 0 0 3652 
16 0 37 490 2732 187 165 31 10 3652 
17 0 2 209 2287 489 496 126 42 3651 
18 0 0 16 1883 689 754 222 87 3651 
19 0 0 0 1300 865 1012 325 150 3652 
20 0 0 0 620 1007 1325 493 207 3652 
21 0 0 0 0 1124 1630 634 263 3651 
22 0 0 0 0 1093 1613 638 307 3651 
23 0 0 0 0 1081 1564 637 369 3651 
Total 139 1511 9361 32825 14129 18545 6899 4232 87641 
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Figure 7. The wind rose at Elmdon (Birmingham Airport) shows the predominance of south-westerly winds, the circle delineating the 14% 
frequencies. The wind roses at Wyton and at Brize Norton have the same characteristic. 





Table V. Wind speed percentage frequencies at Elmdon, 
Brize Norton and Wyton between 08 and 18 UTC. The U 
values are average values in the wind speed ranges. The 
frequencies are similar at the three stations. 





Percentage frequencies 





Wind speed U Elmdon Brize Wyton 
(kn) (ms"') Norton 
Calm 0 5 3 2 
3 1 14 20 11 
4-6 25 25 21 20 
7-10 4.25 32 29 Ke 
11-16 6.75 20 22 27 
17-98 9 4 5 7 





Conditions of high humidity (leading to long, visible 
plumes extending 100 metres or so) are very unusual 
during summer daylight hours; a relative humidity 
greater than 95% during the day occurs less than 4% of 
the time between May and the end of August. 

In winter, however, this increases to around 30% at 
9 a.m. between November and February, and 10% at 
3 p.m. The required cold temperatures would also be 
experienced only in winter for such long plumes. 
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Summary 


The main uses of climate data are listed. Some comments are made on the data requirements for climate 
modelling, detection and diagnostic studies. The discussion concentrates mainly on the data needed to provide 
boundary conditions, initialization and validation for climate models. 


1. Introduction 

Many uncertainties surround our understanding of, 
and ability to model, the climate system. These, and 
other technical limitations, make prediction of the 
nature, magnitude and timing of climate change, 
particularly at scales smaller than global, unreliable at 
present. Increasing concern over the need to detect, 
understand and predict global and regional environ- 
mental changes has highlighted the need to establish 
high-quality observational data sets relating to the 
earth’s atmosphere, oceans, cryosphere and land- 
surface. The basic climate record needs to be better 
established from geological times up to the present and 
decisions must be taken now to ensure more compre- 
hensive global data coverage of the climate system in 
future. 

A prerequisite to designing and implementing a 
Global Climate Observing System (GCOS) is to identify 
the principal uses and users of the data to be collected in 
order to justify such a costly, long-term, internationally 
collaborative operation. 

The major requirements for data in support of climate 
prediction and research, the latter to be interpreted in its 
widest sense, are for: 

(a) monitoring the climate and its variability at 

global and regional scales, thereby enabling quantifi- 

cation of natural climatic fluctuations on a range of 
temporal and spatial scales and the detection of 
climate change, 

(b) attribution of climate change to particular 

causes, as in the search for a fingerprint which will 

enable detection and attribution of climate change 
due to an enhanced greenhouse effect resulting from 
man-made emissions of greenhouse gases, 

(c) detection (and attribution) of the environmental 

impacts of climate change, 

(d) diagnostic studies to elucidate the behaviour of 

the climate system and its component parts, namely, 

atmosphere, oceans, land surface, cryosphere, etc. 





Background paper prepared for the World Climate Research 
Programme ad hoc Workshop on the Planning of the Global 
Climate Observing System, 14-15 January 1991, Winchester, 
England. 


including studies of the mechanisms of natural 
climatic variability, 

(e) development and testing of hypotheses relating to 
local and global climate variations and to the degree 
of predictability of climatic phenomena, 

(f) process studies. Special data are needed in 
support of detailed research studies of a wide variety 
of complex dynamical, physical, chemical and 
biological processes which help govern the state and 
evolution of the climate system. Such specialized 
data sets are likely to need to be highly resolved in 
time and space and therefore gathered for a limited 
period over a restricted area of the globe. High- 
priority process studies include cloud-climate and 
air—sea interactions, 

(g) providing boundary conditions for climate models. 
These include surface properties such as albedo, 
roughness length, vegetation index, soil physical 
parameters, etc. Some such data are also needed to 
identify changes in relevant surface properties, e.g. 
due to deforestation (see also (c) above), 

(h) initialization of climate model integrations, 
especially of the oceanic and cryospheric components, 
and also some of the land-surface hydrological 
characteristics such as the ‘soil moisture content’, 
(i) validation of climate models. A very wide range of 
data, including palaeoclimatic data, are needed to 
assess the performance of the models being used for 
climate simulation and prediction. Model behaviour 
is compared with that of the ‘observed’ climate, often 
leading to further development and improvement of 
the models, and , 
(j) data-assimilation techniques for climate model 
development. Aspects of the climate models can be 
improved through the use of climate observations in 
sophisticated four-dimensional data-assimilation 
schemes, such as those used with state-of-the-art 
operational weather forecasting models. 


Although model and observational studies are 
becoming increasingly interdependent, it is convenient 
for the points to be focused on primarily in this article to 
consider the requirements (a)-(f) above as being, in the 
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main, model-independent, whereas (g)-{j) are directly 
related to the development and use of climate- 
simulation and prediction models. A few comments on 
the data requirements for climate monitoring and 
diagnostic studies are given in section 2. Section 3 
concentrates on the observational needs of climate- 
simulation and climate-prediction models, i.e. items (g)— 
(j) above. 
It is worth stressing that the above needs and uses 
encompass a variety of types of data, for example : 
‘Stand-alone’ instrumental data, including data from 
both in situ and remote-sensing (in particular satellite- 
borne) instruments. Examples of the former are 
surface and upper-air temperatures, rainfall, etc. The 
latter category also includes surface and upper-air 
temperatures, and cloud cover, etc. 


‘Blended’ data, i.e. a data set produced from more 
than one data source; e.g. sea surface temperatures 
derived from both satellite- and ship-borne sensors. 


(Climate studies are by their very nature global in 
character and this necessitates extensive use of space 
observations to provide the necessary coverage and 
continuity. Space observations are, however, in 
many respects difficult to interpret and special 
blending techniques are required in order for them to 
provide reliable information for climate research. In 
particular, they need to be combined with in situ 
observations to resolve ambiguities in their interpreta- 
tion, to provide higher precision, and to enable 
estimates of fields of derived quantities, such as 
surface heat and moisture fluxes.) 


Data assimilated by operational weather forecasting 
models. This is a powerful technique for the optimal 
combination of space-based and in situ observations. 
It provides ‘complete’ global data sets which are 
internally, dynamically and physically consistent, 
e.g. 3-dimensional, global wind fields. However, 
biases may be introduced by changes to the models 
used. 


‘Derived’ data, i.e. data calculated from other directly 
measured or assimilated data, e.g. derived surface 
fluxes of heat, moisture and momentum. Information 
on land-surface characteristics, such as albedo and 
surface roughness, are usually ‘derived’ from a ‘blend’ 
of a variety of data types. 


2. Data requirements for monitoring, 
detection and diagnostic studies 

Studies of climatic processes often require short-lived 
specialized regional-scale data sets and associated field 
programmes. Studies of climate variability, change and 
predictability demand long-time-scale (decades to cent- 
uries or more) data sets which are homogeneous in time, 
i.e. without bias or with a constant bias. Studies of 





global environmental change (including climate change) 
clearly require global data sets, and regional variability, 
change and predictability are being studied increasingly 
in the context of the corresponding global-scale 
characteristics. The climate data sets to be assembled 
must reflect these diverse needs. 

The ‘operational’ and research needs will also often 
lead to a hierarchy of data sets for a given variable, e.g. a 
real-time data set (used for operational weather 
forecasting, and for monthly and seasonal predictions), 
a delayed-mode data set (containing more data than the 
real-time one, and useful for preliminary monitoring 
and research studies), and a slow-mode data set 
(containing all available data with instrumental correct- 
ions, e.g. historical in situ sea surface temperature data 
sets). 

Beside the routine production of analyses for climate 
research using data-assimilation techniques, it is essential 
to continue to process selected in situ and remotely 
sensed observations independently from the integrated 
data-assimilation systems. This is because information 
is needed on the atmosphere and oceans (for example, 
ocean and land surface temperature analyses) which are 
free from the biases that arise from the inclusion of 
(changing) modelling assumptions. For selected variables, 
long time-series of observations need to continue to be 
carefully processed to produce analyses which are as free 
as possible from the influence of changes to the 
observing or processing techniques introduced over the 
period. ‘Pure-data’ archives are necessary for this 
purpose. A particularly important example for validating 
climate models is the need for a reliable, homogeneous 
archive of atmospheric data from radiosondes. 

For monitoring climate, quantifying climate variability 
and detecting climate change, a long historical perspective 
is needed, in contrast to the normal data requirements in 
operational weather forecasting. In this context, it is 
important to realize that many historic data are not yet 
in aform for use on computers. Indeed many of the data 
needed for climate purposes are not expected to be 
processed through data-assimilation schemes in the 
medium-term future. Criteria for the collection and 


assembly of data sets for climate monitoring include : 
Data are required tor different parts ot the climate 


system, particularly for atmosphere and oceans. 
There is a need to establish reliable homogeneous 
quasi-global historic ‘base-line’ data sets. 

Current data need to be homogeneous with historic 
data; this is a difficult task and specialized techniques 
need to be developed to achieve the necessary high 
levels of homogeneity. 

Satellite data need to be blended with more- 
conventional data. 

Near real-time data are required to monitor short- 
period weather variability on time-scales ranging 
from a day to a month. Few daily data are available 
on the global scale and this prevents worldwide 
monitoring of extreme events such as frosts. 
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In the context of short-period weather variability it is 
worth noting the additional data needs for long-range 
(monthly-seasonal) forecasting. This activity requires 
daily, real-time global data sets with special requirements 
for particular regions (e.g. high-resolution climate data 
for the United Kingdom are needed for the long-range 
forecasting activities within the Meteorological Office). 
Also, long-range forecasting based on statistical methods 
needs long, homogeneous historical data sets. 

Some suggested data sets for climate monitoring are 
listed in Table I. Many of the qualifying column entries 
in this table are somewhat tentative. 

Diagnostic studies of data to explore climate mechan- 
isms and processes, especially studies of low-frequency 
weather variability, have some additional and different 
requirements. The overriding need is for 3-dimensional, 
high-resolution, physically and dynamically self-consistent 
global data sets. This is generally only feasible by use of 
data assimilation allied to a global model. Data are 
needed at least twice daily (preferably 6-hourly) to 
avoid biases that can be introduced by diurnal effects. 
Homogeneity is another important requirement, but 
cannot always be achieved. Many of the data sets to be 
analysed are of derived physical and dynamical 
properties of the atmosphere and, where possible, the 
oceans, e.g. diabatic heating rates, eddy fluxes, ageo- 
strophic quantities, etc. 


3. Data needs for climate modelling and 
prediction 

Some comments are offered on the following aspects 
of the data requirements in connection with the 
development and use of global numerical models for 
studying and predicting climate and climate change: 

Boundary conditions. 

Initialization. 

Validation. 

Data-assimilation techniques for model development. 


3.1 Data for model boundary conditions 

This encompasses a wide range of data which are 
prescribed in order to run climate models in various 
configurations (namely, atmosphere-only, ocean-only, 
coupled atmosphere—ocean). They may be held constant 
or updated prescriptively throughout a model integra- 
tion, but they are not determined prognostically or 
diagnostically within the model itself. A range of 
boundary conditions is needed for application at the 
land-surface or air-sea interface, in the upper reaches of 
the atmosphere, throughout the atmosphere, oceans 
and soils, and at the lower ocean boundary. 

The need to prescribe such boundary conditions 
depends to a very large extent on the state of develop- 
ment of the model being used, the configuration it is 
being used in, and what it is being used for. As models of 
the full climate system develop, and as they are used 
in easingly for making predictions of transient global 
and regional climate change, then correspondingly more 


degrees of freedom will be given to such models and the 
requirements for boundary conditions will evolve in 
response. Some quantities acknowledged at present as 
boundary conditions will become full prognostic vari- 
ables in the models, whilst new boundary conditions will 
be needed for new and more complex representations of 
the sub-grid-scale processes (i.e. new parametrizations 
will require new boundary conditions). 

Examples of important boundary conditions requiring 
better global data sets are : 


Atmosphere-only models: 
Sea surface temperatures (SSTs). 
Sea ice extent, concentration, thickness. 
Ozone distribution. 
Land-surface characteristics (albedo, surface rough- 
ness length, vegetation type, etc.). 
Land ice extent and thickness. 
Orographic height (and sub-grid-scale variance). 


Ocean-only models: 

Aside from ocean bottom topography, the prime 
requirement here is for data to provide the surface 
forcing needed for ocean model integrations, namely: 

Surface stress fields. 

Wind mixing. 

Heat fluxes (incoming solar and long-wave radiation, 

sensible and latent heat fluxes) and freshwater 

(precipitation less evaporation) fluxes. 

River runoff. 


Two basic approaches may be made to provide the 
appropriate forcing for the model, though in practice a 
combination of the two is frequently used. 

The fluxes may be specified directly, in which case it is 
necessary to use a ‘flux-correction’ technique whereby 
the net heat flux (Q.) is modified by a feedback term 
which helps to constrain the predicted sea surface 
temperature (7,) in the model to the observed field (7,3) 
depending on the size of the feedback parameter (4), i.e. 
the applied heat flux (Q) is given by: 


Q= Q + k(To — Tp). 


A similar approach may be used for the freshwater 
flux, which is then constrained to the observed salinity 
field. 

Alternatively, the basic meteorological parameters 
necessary to drive the turbulent fluxes can be specified 
and the fluxes derived via appropriate bulk formulae as 
the integration proceeds and with reference to the 
modelled SST (or ice surface temperature) where 
appropriate. In this case the data sets required are: 

Surface wind field. 

Wind mixing. 

Surface air temperature. 

Surface air humidity. 

Surface pressure. 
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Incoming solar and long-wave radiation/ cloudiness. 


Coupled atmosphere—ocean general circulation models 
(AOGCMs): 
Particular boundary data sets required here are: 
Ozone distributions. 
Land-surface characteristics, including land ice extent 
and thickness. 
Topographic data for both land and oceans. 


3.2 Data for model initialization 

These are data needed to set the initial values of 
variables which will then be updated, usually prognostic- 
ally, within the model itself. It is often not understood 
that climate modelling and climate prediction are not 
initial-value problems in the same way that short- to 
long-range forecasting are initial-value problems — at 
least as far as initializing the atmosphere is concerned. 
Initialization of the atmosphere is not a critical factor in 
climate modelling or prediction; it is quite acceptable to 
use model-compatible data from a single time from a 
previous climate model run or a numerical weather 


‘prediction analysis, and quite unnecessary and of no 


additional value to demand real-time global analyses. 

In the current state of development of climate models, 
it may be more critical to initialize aspects of the land- 
surface prognostic variables, if not starting from a 
previously derived model data set. In particular, there is 
modelling evidence of long-lasting effects of the initial 
values used for some of the land-surface hydrological 
characteristics. Snow cover, soil moisture content and 
soil temperatures are examples of variables which may 
need to be initialized by ‘spinning up’ the model through 
at least a seasonal cycle. 

The main current (and probably future) requirements 
for initializing climate models concern the oceans, 
whether in ocean-only models or coupled AOGCMs. 
Global fields of profiles of temperature, salinity and 
currents are required to provide a realistic initializa- 
tion of ocean models. Currently the Levitus global 
temperature and salinity data set is used and the currents 
generated as the model runs. Repeated insertion of 
Levitus data may be useful and valid for high-resolution 
ocean models in this context (such as that used in 
NERC’s FRAM project). Long (centuries or more), 
costly, computer integrations are necessary for the 
models to come into equilibrium with any specified 
forcing. Note that the final state may depart significantly 
from the ‘observed’ state. 

Aspects of the cryosphere also require careful 
initialization in climate models. In particular, as a 
minimum, sea ice extent needs to be initialized. There is 
also need for sea ice concentrations and thicknesses, 
although the latter are particularly difficult to achieve. 

In coupled AOGCMs a major problem is to bring the 
component systems into mutual adjustment with one 
another. The length of model integration necessary to 
achieve sufficient adjustment (and estimate flux correct- 





ions) for a given purpose may vary from decades to 
millenia. The initial data requirements for AOGCMs 
are as above. There is no doubt that initialization of the 
ocean component of the climate system presents by far 
the greatest challenge in this area of climate modelling 
and prediction. Note that careful specification of the 
‘observed’ state is particularly necessary for TOGA- 
related predictions and hindcasts (e.g. the prediction of 
ENSO events). 


3.3 Data for model validation 

The IPCC WGI Report identifies three categories of 
data needed for climate model validation: 

(a) variables important for description of the atmos- 

pheric and oceanic circulation, 

e.g. atmosphere: mean-sea-level pressure; wind and 

temperature profiles; variability as indicated by eddy 

kinetic energy, and 

ocean: surface dynamic height (not often calculated 

in models); temperature, salinity and current structure; 

distribution of tracers; eddy statistics, 

(b) variables critical in defining climate change, 

e.g. atmosphere: surface air temperature; precipitation; 

soil moisture; monthly means plus interannual and 

daily variability, and 

ocean: SST; ocean mixed-layer depth, and 

(c) variables important for climate feedbacks, 

e.g. snow cover; sea ice; clouds and their radiative 

effects. 

Note that the above categories (a)-(c) are not 
mutually exclusive and many types of observations fit 
equally well into more than one. 

Validation of selected regional aspects is also of 
relevance, e.g. occurrence of ENSO events in coupled 
models, and monsoon phenomena. 

There are again particular considerations when 
verifying ocean models. On the seasonal time-scale, data 
are needed to validate simulations of the upper ocean 
mixed layer. Data on the spread of transient tracers 
provide valuable verification of decadal time-scale 
changes in the ocean. On the time-scale of a century 
or more data on deep ocean temperature and salinity 
structures, as well as of other tracers, will be of value. 

Particular requirements for validating results from 
coupled AOGCMs include data on surface fluxes, SST 
and sea ice (extent, concentration and thickness). 


3.4 Additional comments on the needs for 
boundary conditions, initialization and 
validation of climate models 

Table II lists some of the recognized data requirements 
for initializing, validating and providing boundary 
conditions for climate models. 

Data are needed with a range of temporal resolutions: 
monthly means and variances, both on a climatological 
basis and for individual months; daily data are explicitly 
required in some cases as are data on the diurnal cycle of 
some quantities. 
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Table ll. Data required to provide initialization, validation 
and boundary conditions for climate models (tentative) 








Boundary Initialization Validation 

Atmosphere 

T, gq, V (surf & u/a) * (surf) ry 
Cloudiness * . 
Radiative fluxes 

(surface, vertical 

profiles & TOA) * * 
Precipitation * 7 
Surface fluxes * . 
Ozone ¥ * 
Land surface 
Albedo e * 
Snow depth ° ° 
Vegetation e (*) . 
Runoff 2 ® 
Soil moisture . ’ 
Surface temperature ° e 
Topography 7 ‘ 

Roughness e 
Sea ice 

Extent/ concentration * e € 
Thickness * * . 
Roughness ° 
Ocean 
Dynamic height (®) x 
Surface T, S, V ® (SST) * . 
T, S, V structure . ° 
Mixed-layer depth e ° 
Wave characteristics 

(height etc.) . 
Tracers . s 
Water type e 


(plankton, COz, etc.) 





Size of dots signifies probable importance. 
* signifies over ocean only. 


Future model assessments would benefit particularly 
from improved data sets on: 

Precipitation and evaporation rates over the oceans. 

Evapotranspiration, soil moisture and snow depth 

over land. 

Clouds. 

Ocean properties (temperature, salinity, currents, 

ete.). 

Sea ice. 

With regard to the potential of climate data sets 
produced by means of sophisticated data-assimilation 
techniques there is a need for climate modellers to 
explore more thoroughly how ‘operational’ archives can 
be utilized more comprehensively for model validation. 
There is still considerable capital to be extracted from 
existing data sets and scope for carefully considered 
reanalyses of data from the recent past, using a fixed 
model. 

There is increasing need to ensure that data are made 
available in machinable form, compatible with what 
models require. In that context, more-uniform practices 


need to be adopted in the retention of the model data 
which are needed increasingly for model intercomparison 
studies, e.g. snow-cover frequency and depth, extremes 
and means of daily near-surface temperature. 


3.5 Data for assimilation techniques for 
development of climate models 

There are insufficient observations at any one time to 
determine the state of the atmosphere (even less for the 
oceans). Therefore we need to invoke any additional 
information we have, and this is available indirectly as 
the knowledge of the behaviour and structure of the 
atmosphere and, to a lesser extent, of the oceans, which 
we have encapsulated in the formulations of our 
operational weather forecasting and climate models. In 
particular, knowledge of evolution with time is embodied 
in operational forecasting systems and this enables the 
use of data distributed in time. Such models also provide 
a dynamically and physically consistent means of 
representing the atmosphere and oceans, and of deriving 
fluxes and other diagnostic quantities. 

Assimilation is the process of finding the model 
representation which is most consistent with the 
available observations. Given observations distributed 
in time and space, and a forecast model, we can perform 
a 4-dimensional data assimilation. This is normally 
done by adding observations as the model is.integrated 
forward in time. The current model-state summarizes in 
an organized way the information from earlier observ- 
ations and, given new observations, the model state is 
modified to be consistent with these and the earlier 
information. 

It should be emphasized that at any given time the 
model state usually contains more information than can 
be extracted solely from the currently available 
observations. However, only variables which are well 
represented in the weather forecasting model can be 
assimilated sensibly in this way. There is considerable 
need and scope therefore for developing data assimilation 
techniques to process new types of observations within 
the framework of model assimilation in order to 
calculate the misfit between observed (or subsequently 
derived) properties and those produced (or diagnosed) 
from climate model integrations. Such observational 
and model comparisons can be used to validate and 
improve the representation of the physical processes in 
models and, conceivably, when such parametrizations 
have been demonstrated to represent an observed 
property adequately, then the model can be developed 
into a full assimilation using inverse methodology. 
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Awards 


L.G. Groves Memorial Prizes and Awards 
for 1989 

At ajoint ceremony at RAF Odiham on 10 December 
1990 The Chief of Air Staff, Air Chief Marshal Sir Peter 
Harding, GCB, ADC, D.Sc., F.R.Ae.S., CBIM, presented 
the Wilkinson Battle of Britain Memorial Sword to the 
Support Helicopter Force and the L.G. Groves Prizes 
and Awards to the prize-winners. Air Vice-Marshal 
D.O. Crwys-Williams, RAF, Director General of 
Personnel Services, Air Officer Commanding No.1 
Group, read the citations and announced the prize- 
winners. 


Meteorology Prize — A.C. Lorenc 





The citation for this award was: 

“Mr Lorenc is one of the world’s leading experts on 
research into effective methods of processing meteor- 
ological observations automatically to provide initial 
conditions for numerical weather forecasts. A hallmark 
of his work is its sound theoretical basis and its clear 
practical applicability. 

Recently he showed that quality-control procedures 
could be built upon well-found statistical principles, and 
he demonstrated that automatic techniques can supplant 
time-consuming human scrutiny of data for errors. 

Another of his innovations is the scheme currently 
used to provide operational analyses. One of its features 
is the ability to assimilate asynoptic observations (such 
as aircraft data) at the actual time that the measurements 
were made rather than (as was previously done) at the 
nearest 6-hour time slot; such a procedure makes better 
use of the data especially in rapidly developing synoptic 
situations. Its operational implementation in late 1988 
has lead to improvements in the accuracy of wind 
forecasts for aviation. 

More recently his attention has turned to the effective 
use of remotely sensed data. He is Principal Investigator 
for the American Earth Observation Satellite and for 
the European ERS-I satellite. His advocacy of the use of 





numerical model forecasts, rather than climatology, asa 
first guess for satellite-derived temperatures has led to 
improved forecasts, and has strongly influenced US 
meteorologists to adopt a similar approach.’ 


Meteorological Observation Award — 
W.D.N. Jackson 





The citation for this award was: 

“Mr Jackson has been at the Meteorological Research 
Flight for 34 years. During this time he has been in 
charge of the team of about 20 people responsible for the 
aircraft instrumentation and the airborne and ground- 
based computer systems used to analyse the resulting 
data. He has caused significant improvements and 
advancements in all these areas and during 1990 a new 
display and analysis system was completed. Large parts 
of these systems were designed and implemented by 
Mr Jackson himself. The end result is much better 
aircraft utilization due to both the new real-time 
displays of data on the aircraft and the ability to analyse 
flight data in detail overnight so as to identify problems 
and issues prior to subsequent flights. In addition to 
these technical advances Mr Jackson has created an 
atmosphere in which the staff recognize the importance 
of aircraft availability and have a framework within 
which to strive to maximize it. He has made a most 
significant contribution to the observing programme of 
the Meteorological Research Flight.’ 


Other prizewinners 

The Air Safety Prize was presented to Fit. Lt. 
R.L. Foulkes B.Sc., RAF, but the Ground Safety 
Award was not presented as the prize-winner, Sqn. Ldr. 
M.A. Calame, RAF, was unable to attend. 
Correction 
Meteorological Magazine, April 1991, p. 73, Review of 
Global Air Pollution. 
The specific reference to Professor Scorer’s book was 
accidentally transposed to the second paragraph, but 
should have been included at the end of the first 
paragraph, as an alternative book to the one under 
review. Apologies are offered for any implied mis- 
representation. 
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Satellite photographs — 2 February 1991 at 1533 UTC 


The NOAA-11 pictures exhibit features which have 
come to be recognized as hallmarks of the onset of 
explosive cyclogenesis. 

The images display several prominent features as 
highlighted on the schematic diagram in Fig. 1: 

(a) Cloudy band, F, which is symptomatic of a warm 

conveyor belt flow associated with a well-established 

cold front. 

(b) The ‘cloud head’ or ‘baroclinic leaf’, C, which is 

particularly striking. It is interpreted as a region of 

strong slantwise ascent occurring rearwards of a 

newly developing front RR. The corresponding 

infra-red image reveals a transition from warm cloud 
tops at the eastern edge of this cloud shield to coldest 
cloud tops near its sharply defined western boundary 

— a feature typical of slantwise ascending motion. 

(c) A narrow cloud-free region separates cloud 

masses C and F — this is particularly prominent on 

infra-red imagery. Such ‘dry intrusions’ or ‘dry 
wedges’ (see Monk and Bader*) are often a feature of 
rapidly deepening depressions. 

Further to the west, streets of convective cloud, P, 
stream eastwards around an upper trough, the accompany- 
ing cold advection reinforcing the thermal gradient 
along the new frontal zone. Lines of convective cloud 
can be seen undercutting the overhang of cirrostratus. 


The southern, cyclonically curved, tail of the cloud 
head consists only of low cloud, at the leading edge of 
which is a well-marked narrow ‘rope cloud’ (RR). This 
probably represents line convection at the surface cold 
front which has been overridden by the intrusion of dry 
air aloft. 

The transverse striations in the cloud head have been 
observed previously, notably in the Great Storm of 
October 1987 (see ShuttsT), but are not fully understood. 
The latitudinal banding became less pronounced in later 
satellite passes, obscured by thickening cirrus layers 
generated by continued slantwise ascent. 

The developing low, L, was estimated to have been 
centred near the northern end of the rope cloud, with a 
central pressure of about 990 mb. The midday operational 
analysis of the fine-mesh model had identified only an 
open wave of about 1000 mb, but 24 hours later the low 
had deepened to below 935 mb as it crossed Iceland, 
causing extensive structural damage. Many Icelandic 
stations observed hurricane force winds: Vestmannaeyjar 
on the south coast reported a mean wind of 103 kn at 
1200 UTC on 3 February. 

Although the numerical weather prediction guidance 
based on data at 1200 UTC on 2 February failed to 
capture the intense degree of cyclogenesis, the imagery 
was used as a basis for adjusting later model analyses to 





Figure 1. 


Schematic diagram of the features described in the text. 





* Monk, G.A. and Bader, M.J.; Satellite images showing the develop- 
ment of the storm of 15-16 October 1987. Weather, 43, 1988, 
130-135. 


+ Shutts, G.J.; Dynamical aspects of the October storm 1987: A study 
of asuccessful fine-mesh simulation. Q J R Meteorol Soc, 116, 1990, 
1315-1347. 
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Figure 2. Infra-red image of the situation shown schematically in Fig. - Photograph by courtesy of University of Dundee 





Figure 3. Visual image of the situation shown schematically in Fig. 1. Photograph by courtesy of University of Dundee 
cater for the suspected existence of a much deeper 
depression. This successfully culminated in a much 
improved analysis, and a subsequent operational run of 


the model subsequently deepened the depression to 
946 mb accompanied by an extremely vigorous circulation. 


J. Norris and M.V.Young 
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